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GALACTIC HYDRAULIC DROP? 



STEEN H. HANSEN 
University of Zurich, Winterthurerstras.se 190, CH-8057 Zurich, Switzerland 



Galaxies and galaxy clusters are observed to have a rather non-trivial radial behaviour. The 
observations show that the radial profiles change from one power-law profile near the centre to 
another power-law profile in the outer region. We present a simple explanation for this complex 
behaviour by finding the analytical solutions to the governing hydrodynamic equations. We 
see that the origin of this complexity is the collisional nature of the baryonic plasma, possibly 
related to a turbulence-enhanced viscosity. 



1 Introduction 

Large gaseous baryonic structures such as galaxies and galaxy clusters have been known and 
observed for many years. A characteristic behaviour is that the radial profiles, e.g. of surface 
brightness or electron density, often have the complex behaviour that they follow one power- 
slope, a, in the inner part, and another power-slope, j3, in the outer part beyond a characteristic 
radius, ro, 

(i) 



A 



(*) ( 



1 + 



ro 



and this transition is frequently observed to be rather sharp. This behaviour is often fit by 
observers by simple phenomenological profiles like eq. (J2) which is composed of just such two 
power-laws. These observations include spiral galaxies (e.g. using WFPC2 data 1 z ) and clusters 
of galaxies (e.g. using Chandra data^). However, there is little (if any) theoretical guidance for 
the use of such profiles. Here we attempt a derivation of this complex behaviour. 

Traditionally one would expect that two different power-slopes must be related to different 
physics, in particular if the transition is sharp. Surprisingly enough, for the radial density profile 
this does not have to be the case. We will show that the governing equations have exactly two 
solutions, which imply that the inner and outer density profiles may choose different solutions, 
and hence quite generally will be different (see ref. [4] for details). 




Figure 1: The circular hydraulic jump, as observed in any kitchen sink. The water chooses one solution inside the 
jump and another solution outside the jump. We show that galaxy clusters behave in a similar manner, and that 
the density profile (or surface brightness) therefore is expected to exhibit a break at some characteristic radius, 

providing us with a galactic hydraulic drop. 

In the field of hydrodynamics cases are known where a given set of equations have two 
solutions, and that Nature chooses to use both solutions simultaneously. One well-known example 
is the hydraulic jump, which is a centimetre large ring, observed in any kitchen sink when the 
water flows out radially after hitting the sink (see Figure 1). The water in the in ner f ew 
centimetre follows one solution, and outside the jump the water follows another solutioJ3^. It 
turns out that in a similar manner the density profile in the inner part of e.g. a relaxed galaxy 
cluster chooses one solution, whereas the outer part of the same cluster chooses another solution. 

2 Solving the Navier-Stokes equations analytically 

The behaviour of any collisional gas or fluid is fully determined by the Navier-Stokes (N-S) 
equations. Baryons often have sufficient collisions to be described by the N-S equations, e.g. in 
a typical intra-cluster gas the equilibration timescale is about 10 years, with mean free path of 
tens of kpc compared to radii of few Mpc. 

In this proceeding we will for simplicity consider a stable spherical cluster of galaxies, where 
the system has picked out an orientation in space, such that all the gas is moving only in the 
0-direction. Thus we have v r = vj, = 0. Here we use notation where r is the radial coordinate, 
O is the angle in the xy-plane, and (ft is the angle from the z-axis. One must keep in mind 
that by considering the N-S equations we are taking a fluid approach which implies that we 
are following a fluid element, and this basically corresponds to averaging over all the particles 
moving through that fluid element. For the ©-velocity we consider the simple form 



where a is the constant to be determined first, v a and r a are unknown constants, with the 
physical interpretation that r a is a characteristic transition radius, and v a is the velocity of the 
fluid element at that radius. 

The first N-S equation. The first N-S equation becomes very simple with the assumed 
form of the velocities 




(2) 



= i/ V 2 u 



(3) 




where V 2 is the scalar Laplacian, and v is the kinematic viscosity. For now all that matters is 
the existence of a non-zero viscosity, so the absolute magnitude (and even radial dependence) is 
not important for the results. Certainly baryons have non-zero viscosity, however the viscosity 
could also arise from turbulence, in which case it could be very larg 

^turb ~ lAv ' wnere tne 

dimension I is the size of the turbulent eddies, and Av is the velocity dispersion. Numerical 
simulations^ have shown that such turbulence indeed exists in galaxy clusters with Av ~ 
300 — 600 km/sec, and I ~ 100 — 500 kpc, leading to a very large (turbulence enhanced) viscosity. 
When we use the form for vq in eq. (J2J), then eq. © has exactly 2 solutions 

a = 1 or - 2 . (4) 

Thus, looking at eq. (|2*|). it is clear that the solution with a = —2 is divergent for r —* 0, and 
we will therefore refer to this solution as the 'outer solution', and similarly, the solution with 
a = +1 is inconsistent for large radii, and we will refer to this as the 'inner solution'. Hence the 
general flow pattern changes from the inner to the outer region, and one may therefore expect 
to find different density profiles in the central and outer regions. 

The second N-S equation. We can now use the next N-S equation to extract the asymp- 
totic radial density profiles. Also this equation is very simple 

v% _ ldP M(r)G 

— ^ 5 — ) (5) 

r p or r z 

where p is the radially dependent density, P is the pressure, G is the gravitational constant, 
and M(r) is the mass within the radius r. We assume that the pressure and density are related 
through P = P a ifij pa) 1 , where P a and p a are the unknown pressure and density at r a . We 
take a monatomic gas with 7 = 5/3. Let us consider densities of the form 

p{r)=pa (£)"' (6) 

such that the parameter (3 determines the density profile. It is worth emphasizing that it is 
exactly this (3 which we are trying to find. 

The last (gravitational) term including M{r) depends on the given system we are considering. 
If the mass is dominated by a point gravitational source (e.g. a central black hole (BH)), then 
it goes like M(r)G/r 2 ~ r~ 2 . If the mass is dominated by the matter density, then it goes like 
M(r) ~ J p(r)dV, with dV the volume element. For spherical solutions this gravitational term 
thus goes like r^ +1 with (3 from eq. ©• 

It is now straight forward to solve eq. (JHJ) in the inner and outer regions. The a = 1 is 
the inner solution, for which we find (3 = —6 when the baryons dominate the mass. With 
BH dominance one has (3 = —3/2, and if dominated by another spherical distribution (which 
probably should arise from dark matter (DM)) with profile (3 S , then we find [3 = 3/2 (f3 s + 2). 
For a different polytropic index, 1 < 7 < 5/3, the coefficient changes from —3/2 to (7 — 1) ■ 
Thus, if the DM has a slope of /3 S = —2 (as expected from adiabatic contraction), then the 
baryons develop a core, (3 = 0. 

In the outer region (a = —2) we find f3 = —6. If a BH dominates then (3 = —3/2, and 
again if a another spherical distribution (DM) dominates then 2>/2{(3 s + 2). E.g. a DM slope 
of —3 leads to (3 = —3/2. One should, however, keep in mind that in the outer region there 
are possibly not sufficient collisions to assure the validity of the N-S equations, so those simple 
solutions should not be trusted too much. 



3 How to test these findings 



Our main finding is that a tran sition from the inner to the outer region generally exists. This is 
observationally well established^^El. \y e have thereby provided theoretical support for the use 



of phenomenological profiles like eq. lfl|). One can now take a further step and test the actual 
numbers we find in our simplified treatment. 

The findings for BH domination always give (5 = —3/2, which is just what numerical simula- 
tions find^l The actual density profiles can be observed in different ways. X-ray observations 
of the luminosity in various bands give the electron density of the plasma as a function of radius. 
E.g. for the relaxed cluster A2029 the outer baryonic profile of (5 = —1.62 was found^. In the 
future the Sunvacv-Zcldovich effect will directly provide a measure of the radial electron den- 
sity 1 ' v \ which will probe large cluster radii since the SZ effect is proportional to n e whereas 
X-ray observations are proportional to r? e . Surface brightness from radio observations of HI 
and molecular gas can in principle determine the baryon profile^l For details and discussion 
see ref. [4]. We are looking forward to doing a more detailed analysis. 



4 Conclusions 



We have presented an explanation for the origin of the complex radial structure of galaxy clusters. 
Specifically, we have shown that the density profiles generally are expected to make a transition 
from one power-slope in the inner to another power-slope in the outer region 

Ps«( r ) = r/ 9i(x + r)^ 2 ' ^ 

The physical origin of this complexity is the collisional nature of the baryonic plasma, and we 
speculate that it may be related to a turbulence-enhanced viscosity. 
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